AD-AT28  397 
UNCLASSIFIED 


FORTRAN  IMPLEMENTATION  OF  COMPLEX  LEAST  SQUARES 
ADAPTIVE  LATTICE  ALGORITH. . (U)  PENNSYLVANIA  STATE  UNIV 
UNIVERSITY  PARK  APPLIED  RESEARCH  LAB.. 

B  A  COOPER  ET  AL .  12  NOV  82  F/G  12/1 


>6-83 


OTIC  FILE  COPY  AD  4  1 2 


#3 


(£) 


FORTRAN  IMPLEMENTATION  OF  COMPLEX  LEAST  SQUARES 
ADAPTIVE  LATTICE  ALGORITHM  AS  AN  ALL  ZERO  INVERSE 
FILTER 


B.  A.  Cooper  and  L.  H.  Sibul 


Technical  Memorandum 
File  No.  TM  82-219 
November  12,  1982 
Contract  No.  N00021-79-C-6043 


Copy  No.  30 


Copy  available  to  DTIC  does  not 
permit  fully  legible  reproduction 


The  Pennsylvania  State  University 
Intercollege  Research  Programs  and  Facilities 
APPLIED  RESEARCH  LABORATORY 
Post  Office  Box  30 
State  College,  PA  16801 


Approved  for  Public  Release 
Distribution  Unlimited 


NAVY  DEPARTMENT 


NAVAL  SEA  SYSTEMS  COMMAND 


DTIC 

^ELECTEJ* 

MAY  2  3  19831  ■ 

E 


°s  23  0Q8 


I 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  rage  I'W'han  Oma  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

i  report  number 

TM  82-219 

WBBmM 

3  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

FORTRAN  IMPLEMENTATION  OF  COMPLEX  LEAST  SQUARES 
ADAPTIVE  LATTICE  ALGORITHM  AS  AN  ALL  ZERO  INVERSE 
FILTER 

5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

Interim 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  authors; 

B.  A.  Cooper  and  L.  H.  Sibul 

8.  CONTRACT  OR  GRANT  NUMBERfsJ 

N00024-7  9-C-6043 

9.  performing  organization  name  and  aooress 

Applied  Research  Laboratory 

P.  0.  Box  30 

State  College,  PA  16801 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  &  WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

D.  Houser,  Code  63R-14 

Naval  Sea  Systems  Command 

Department  of  the  Navy,  Washington,  DC  20362 

12.  REPORT  DATE 

November  12,  1982 

13.  NUMBER  OF  PAGES 

51 

l«  MONITORING  AGENCY  NAME  6  A00RESS(7f  dlllerent  Irom  Controlling  Oltice) 

15.  SECURITY  CLASS,  (of  this  report) 

Unclassif ied 

IS*.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (o  I  this  Report) 


Approved  for  public  release,  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  the  abetrect  entered  In  Block  20,  II  dlllerent  from  Report ) 


18  SUPPLEMENTARY  notes 


I 


19.  KEY  WORDS  (Continue  on  reverse  aide  if  necessary  and  identify  by  block  number) 


fortran,  adaptive,  algorithm,  speech,  signal  processing 


20.  ABSTRACT  (Continue  on  reverse  side  it  necessary  and  identity  by  block  number) 


x  A  complex  least  squares  adaptive  lattice  algorithm  (1)  has  been  imple¬ 
mented  on  the  ARL  VAX/VMS  computer  as  an  inverse  all  zero  filter.  A  FORTRAN 
program  package  is  presented  which  also  includes  a  recursive  realization  of  a 
pole-zero  prefilter  and  an  input  signal  generator.  A  series  of  examples  is 
taken  from  the  speech  field,  using  the  inverse  filter  to  deconvolve  speech 
signals  produced  by  vocal  tract  models  of  varying  complexity.  The  input  to 
the  prefilter  is  shown  to  be  accurately  recovered  as  the  output  of  the  inverse 
filter.  . 


DD  I  JAN  7J  1473  EDITION  OF  I  NOV  «S  1$  OBSOLETE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  this  PAGE  ftHi»n  Date  Entered 


id 


, '■’cession  For 

'•'TtT  gra&i  ^ 

me  tab  □ 

Unannounced  Q 

Justification - 

I  By__ - - - 

Distribution/ 

Availability  Codes 
Avail  and/or 
Dist  I  Special 


Subj  ect : 


Re  ferences : 


ABSTRACT 


FORTRAN  Implementation  of  Complex  Least  Squares  Adaptive 
Lattice  Algorithm  as  an  All  Zero  Inverse  Filter 

See  Pages  14  and  15. 


A  complex  least  squares  adaptive  lattice  algorithm  (1)  has  been 
implemented  on  the  ARL  VAX/VMS  computer  as  an  inverse  all  zero  filter. 
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series  of  examples  is  taken  from  the  speech  field,  using  the  inverse 
filter  to  deconvolve  speech  signals  produced  by  vocal  tract  models  of 
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INTRODUCTION 

A  complex  least  squares  lattice  adaptive  algorithm  (1)  has  been 
implemented  on  the  ARL  VAX/VMS  computer  as  an  inverse  all  zero  filter. 

An  interactive  FORTRAN  program  package  is  presented  which  also  includes 
a  recursive  realization  of  a  pole-zero  prefilter  and  an  input  signal 
generator.  No  knowledge  of  FORTRAN  is  required  to  use  the  package.  The 
input  signal  characteristics,  filter  coefficients  and  critical  lattice 
parameters  are  all  user  specified. 

Three  classes  of  adaptive  algorithms  are  treated  in  the  literature: 
(Widrow's)  LMS ,  the  gradient  lattice,  and  the  least  squares  lattice.  Of 
the  three,  the  least  squares  adaptive  lattice  algorithm  has  been  shown 
to  have  the  fastest  convergence  rate  (2,3).  For  this  reason,  it  has 
been  used  successfully  in  a  wide  range  of  applications;  i.e.,  data  decon¬ 
volution,  system  identification,  speech  processing,  and  spectral  whitening. 
The  program  package  described  in  this  memorandum  was  developed  with  these 
applications  in  mind. 

BACKGROUND 

Consider  a  filter  with  transfer  function 


H  (z) 


Z(z)  _  B(z) 
V(z)  A(z) 


(1) 


where 

NPOLii  NPOLE 

A(z)  =  1  +  Z  A  z  =  l  A  z  ,  A  =  1 
k=l  k=0  k  0 
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and 


NZERO  NZKRO 

B(z)  =  1  +  Z  B  z  =  Z  B.z  ’  B0 
k=l  k=0 


then  Z(z)A(z)  =  B(z)V(z)  and,  after  taking  inverse  z  transforms  of  both 
sides,  the  filter  can  be  realized  in  the  time  domain  by  the  following 
recursion: 


NPOLE 


NZERO 


z(n)  +  Z  A,  z(n-k)  =  vfn)  +  Z  B,  v(:i~k) 
k=l  k  k=l  k 


NPOLE 

z(n)  =  -  I 
k=  1 


z(n-k)  +  v(n)  + 


NZKRO 

r 

k=  1 


vfn-K) 


The  roots  of  A(z)  are  the  poles  of  the  svstem;  the  roots  of  B(z)  are  the 
zeroes  of  the  system.  For  stability,  these  roots  must  lie  outside  the 
unit  circle  because  they  are  roots  of  polynomials  in  z  F . 

The  resonances  and  antiresonances  of  the  human  vocal  tract  (for  a 
particular  steady  state  vowel  sound)  can  be  modeled  by  poles  and  zeroes, 
respectively,  in  the  z  plane.  It  is  a  common  practice  in  the  speech  field 
to  model  these  formants  with,  an  all  pole  model  rather  than  with  poles  and 
zeroes;  adding  a  sufficient  number  of  poles  has  been  shown  to  closely 
approximate  the  effect  of  the  zeroes.  The  roots  of  the  denominator  poly¬ 
nomial  are  related  to  the  vocal  tract  formants  by  the  following  express¬ 
ion  (5)  : 

*  -°kT  ±j2TTFkl 

V  zk  z2k’  z2k-l  e  e 
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where  F,  is  the  formant  frequenc/  and  2 O,  is  the  bandwidth  of  the  kt!l 
k  k 

formant.  There  will  be  a  pair  of  complex  conjugate  poles  for  each  formant 
frequency;  the  order  of  the  system  will  be  twice  the  number  of  formants. 
Then 


H(z)  = 


NPOLE 

7T  (1+z  z  ) 
k=l 


For  the  all-pole  model,  the  time  domain  recursion  becomes 


NPOLE 

z(n)  =  -I  A  z(n-k)  +  v(n). 
k=l  k 


An  inverse  filter  for  the  system  of  Equation  1  would  have  the  trans¬ 
fer  function 


H(z) 


A(z) 
B(z)  • 


For  the  speech  case  (all-pole  prefilter),  the  inverse  filter  would  be 
all-zero : 


H 

pre 


(z) 


1 

A(z) 


H. 

inv 


(z) 


1 


H 

pre 


(z) 


A(z)  . 


Obviously,  if  the  systems  are  cascaded,  the  input  to  the  prefilter  will 
be  recovered  as  the  output  of  the  inverse  filter. 

A  complex  adaptive  lattice  structure  has  been  used  to  develop  the 
inverse  filter.  The  actual  lattice  algorithm  is  taken  from  (1).  Users 
are  referred  to  (2)  and  (3)  and  the  references  cited  therein  for  the 
derivations  of  the  complex  least  squares  adaptive  lattice  algorithm. 


-7- 


November  12,  1982 
BAC-I.HS  :cac 


PROGRAM  PACKAGE 

The  lattice  program  package  consists  of  the  following  interactive 
FORTRAN  programs: 

GEN. FOR  input  generator 

FILTER. FOR  prefilter  and  realization  of 

the  lattice  filter 

LSL.FOR  complex  adaptive  lattice  processor 

the  following  subroutines: 


SYNTH 

synthesizes  complex  polynomial 
from  formants  and  bandwidths 

SYNTH1 

synthesizes  complex  polynomial 
from  its  roots 

ZCPOLY 

solves  polynomial  with  complex 
coef  f icicnts 

following  functions: 

MTH$ RANDOM 

(VAX/ VMS  library)  uniform  random 
number  generator 

AMAXLOC 

finds  location  of  highest  value 
in  an  array 

ARRM.'VX 

finds  highest  value  in  an  array 

The  execution  of  the  programs  may  be  controlled  by  either  one  of 
two  command  procedures:  LATTICE.COM  when  using  a  VT100  terminal,  or 
LATTEK.COM  when  using  the  Tektronix  terminal  in  Room  358,  ASB.  These 
command  procedures  are  given  in  Appendix  A.  They  may  be  invoked  by 
typing  ^LATTICE  and  FLATTER,  respectively.  The  following  command  should 
be  placed  into  either  the  command  procedure  or  the  user's  LOGIN.COM  file 
to  allow  access  to  the  IMSL  package.  If  LNK$LIBRARY  has  already  been 
assigned,  use  LNKSLIBRARYl  or  the  next  highest  unassigned  library. 

$ASSICN  DRAO: [SYSTEM] IMSLIBS .OLB  LNK$LIBRARY 


» *-  .■ 
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Appendix  B  presents  a  sample  terminal  session  tor  the  case  illustrated 
in  Figure  3.  For  the  user  interested  in  modifying  the  programs,  the  rest 
of  this  memorandum  deals  with  the  structure  of  the  individual  programs. 

Figure  7  illustrates  the  calling  hierarchy  of  the  lattice  package. 

The  inputs  to  program  LSL  are  Gaussian  white  noise  and  the  output  time 
series  z(n)  from  program  FILTER. FOR.  The  outputs  from  LSL  are  the  set  of 
inverse  filter  coefficients  (A^;  k=l,...  NZERO)  to  which  the  lattice  has 
adapted  after  the  requested  number  of  iterations,  and  the  time  series 
output  ep(n)  which  is  equivalent  to  the  result  of  passing  the  inp  time 

series  z(n)  through  the  transfer  function  H(z)  =  A(z) .  The  outpi  'liter 

coefficients  are  fed  back  into  program  FILTF.R  to  be  realized  reci  vely 
with  input  z(n).  Thus,  since  program  FILTER  operates  for  any  H(z  .c 
is  both  a  prefilter  and  the  realization  of  the  lattice  filter. 

A  sample  case,  then,  may  be  illustrated  by  a  series  of  five  plots 
(Figures  3-6).  The  relationships  between  the  plots  are  shown  by  the 
block  diagram  in  Figure  2: 

(a)  input  time  series 

(b)  prefilter  transfer  function 

(c)  output  of  the  prefilter 

(d)  lattice  (inverse)  filter  transfer  function 
(realized  by  program  FILTER) 

(e)  output  of  the  lattice  filter 

So  that  the  filter  outputs  may  be  most  easily  evaluated,  the  filter  outputs 
(c)  and  (e)  for  deterministic  prefilter  inputs  are  time  series  plots;  for 
a  stochastic  input,  the  outputs  (c)  and  (e)  are  power  spectra.  The  time 
series  plots  (c)  and  (e)  (from  above)  for  the  deterministic  case  are 
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presented  in  Figures  3-6  ior  one  period  of  the  waveform  after  the  time 
series  has  stabilized  (or  adapted).  Again,  the  output  of  trie  Lattice 
filter  should  approximate  the  input  to  the  prefilter  if  tire  lattice  is  a 
true  inverse  filter.  It  is  obvious  that  some  measure  of  goodness  of  the 
inverse  filter  should  be  applied  to  the  output  of  the  lattice  filter. 

At  present,  only  a  graphical  comparison  has  been  attempted.  However, 
dray  and  Market  (7)  suggest  a.  spectral  flatness  measure  which  will  help 
to  smooth  the  discontinuities  in  the  power  spectrum  caused  by  the  occur¬ 
rence  of  the  inverse  filter  zeroes  being  offset  by  a  few  Hz.  from  the 
pretilter  poles.  Also,  a  mean  square  error  criterion  might  be  applied  to 
the  difference  between  the  deterministic  time  series  output  (error  signal) 
and  tiie  input  time  series.  The  transfer  function  of  the  lattice  filter 
should  be  the  inverse  of  the  prefilter  transfer  function.  Thus,  where 
the  prefilter  has  poles  (peaks)  the  inverse  filter  will  have  zeroes  (dips) 

GEN ,  FPK 

Program  GEN. FOR  generates  both  Gaussian  white  noise  ana  an  impulse 
sequence.  A  more  accurate  approximation  to  the  glottal  wave  is  given  in 
(7).  However,  since  the  thrust  of  this  work  is  aimed  at  demonstrating 
that  the  lattice  filter  is  indeed  an  accurate  inverse  filter,  the  two 
cases  used  so  far  are  sufficient  models  of  voiced  and  unvoiced  steady 
state  vocal  tract  excitations. 

Program  GEN. FOR  calls  the  uniform  rand or  number  generator,  MTHSRANDOM 
from  the  VAX/ VMS  library.  The  Cans  si  an  noise  is  then  developed  following 
a  technique  from  (9).  The  output  time  series  are  in  complex  format  but 
the  imaginary  parts  are  zero.  Ttie  sampling  frequency  is  10  kHz.  For 
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the  complex  impulse  input,  the  frequence  is  variable;  the  examples  shown 
in  Figures  J  and  5  use  FO  =  125  Hz,  a  reasonable  fundamental  frequency 
for  an  adult  male  voice.  The  interactive  user  inputs  to  program  GEN 
are  FO,  fundamental  frequency;  M,  number  of  desired  samples;  and  ~,  the 
standard  deviation  of  the  noise.  The  time  series  are  each  stored  in  an 
output  disk  file. 

FILTER. FOR 

The  user  selects  the  desired  type  of  input  and  the  source  of  the 
filter  coefficients.  Several  options  exist  concerning  specification  of 
coefficients.  The  formant  frequencies  and  their  bandwidths  may  be  spec¬ 
ified  and  the  coefficients  calculated  from  subroutine  SYNTH  and  then 
stored  on  disk.  The  stored  coefficients  may  be  used  with  another  input 
opt  ion, read ing  in  coefficients  from  disk.  The  coefficients  nay  be  typed 
in  directly  from  the  keyboard.  For  the  realization  of  the  lattice  (inverse) 
filter,  the  forward  predictor  (outputs  from  program  LSL)  coefficients 
may  be  read  in  from  disk  as  the  lattice  filter  coefficients.  The  other 
interactive  user  inputs  to  FILTER. FOR  are  M,  desired  number  of  time 
samples;  NPOLE,  number  of  poles;  and  N'ZERO,  the  number  of  zeroes.  It  is 
possible  to  specify  a  prefilter  transfer  function  which  lias  both  poles 
and  zeroes  (the  inverse  filter  would  still  be  all  zero).  The  modeling 
of  a  pole  zero  process  by  an  all-pole  process  is  discussed  by  Friedlander 
and  Maitra  (6).  The  examples  presented  here  involve  prefilters  containing 
only  poles  for  the  purpose  of  illustrating  the  accuracy  of  the  all-zero 
inverse  filter. 


»*  ♦  —  —  *  T  *  • 
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C.  =  10  LOG 


where  a  =  0  and  a“  =  a“  =  a“  +  a-  . 

Z  X  V 


•d  by 

the 

f  o  1  low  i  iv, 

M 

9 

I 

z~  (n) 

> 

n=  1 

out 

M 

2 

V 

v  f  n ) 

in 

n=l 

') 

The  filter  transfer  function  is  realized  from  the  filter  coefficients 
by  the  following  formula: 


20  LOG  (HU)‘  =  20  LOG  --J— !  =  20  LOG 

A ( z )'  I  A  t : 


1  .  _1 _  ‘ 

AU)  A*(z)  j 


tor  the  prefilter  and 


20  LOG  |  A(z)  |  =  20  LOG  (A  (.  z)  •  A"  (z  ) )  ^ 
for  the  lattice  filter,  where 

NPOLE  .  . 

A(z)  =  1  +  Z  A.  . 

k=l  k 

The  frequency  range  of  the  plots  is  0-5000  Hz. 

The  output  power  spectrum  is  computed  (4)  for  the  stochastic  case 
from  10  LOG  P(oj)  where 

p(-)-p.„(— )-[.mI2-p„(.)-a(„.1a,(z)  • 

2  2 

For  the  prefilter,  P.  (uj)  =  a.  =  o  .  For  the  lattice  filter,  P.  (jj) 
m  in  v  in 

2 

will  be  P  (^)  from  the  prefilter,  and  P(w)  =  P.  (to)  •  |H(.o)|  = 

out  in 

P i n ( to ^  •  A(z)  •  A"(z).  Again,  the  frequency  range  of  the  plots  is  0-5000  Hz 


a* '<>■?. 
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LSI,  ■  FOR 

Figure  1  shows  the  structure  of  the  least  squares  complex  adaptive 
lattice.  Table  I  lists  the  lattice  parameters  and  their  relationships  to 
the  text  equations  in  (3)  and  to  the  FORTRAN  code  in  (1).  The  algorithm 
from  (1)  has  been  modified  (from  RATFOR)  for  use  by  the  ARL  VAX/VMS 
computer.  Users  are  referred  to  (2)  for  a  detailed  explanation  of  filter 
parameters.  The  interactive  user  inputs  to  program  LSL  are  SNR,  desired 
signal  to  noise  ratio  in  dB;  ALSL,  value  of  orc^er  ° f  the  filter; 

and  LIM,  the  desired  number  of  iterations.  After  LIM  number  of  itera¬ 
tions,  the  forward  predictor  coefficients  become  the  filter  coefficients 
of  the  inverse  filter.  The  convergence  properties  of  this  lattice 
structure  have  been  studied  by  Hodgkiss  (2,3)  and  are  not  considered 
here.  The  maximum  allowable  number  of  iterations  is  set  at  5000  which 
allows  the  lattice  more  than  enough  time  to  adapt. 

Subroutine  ZCPOLY  solves  the  filter  polynomial  to  make  sure  the 
roots  are  outside  the  unit  circle.  If  a  root  is  inside  the  unit  circle, 
it  is  inverted  and  a  new  polynomial  is  synthesized  with  subroutine  SYNTH1. 
The  coefficients  are  then  stored  on  disk  for  use  by  program  FILTER. 

The  forward  and  backward  reflection  (PARCOR)  coefficients  are  also 
stored  on  disk  for  possible  future  use.  These  parameters  are  related  to 
the  reflection  coefficients  of  an  acoustic  tube  model  of  the  vocal  tract. 
The  next  step  in  this  research  effort  will  involve  making  use  of  these 
coefficients  to  directly  identify  the  prefilter  transfer  function. 

Table  II  lists  the  prefilter  coefficients,  final  lattice  coefficients, 
and  both  sets  of  final  PARCOR  coefficients  for  the  four  sample  cases 
presented  in  Figures  3-6. 
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The  desired  signal  to  noise  ratio  is  achieved  by  scaling  down  the 
added  noise  (of  user  specified  standard  deviation  from  GEN'. FOR)  rather 
than  scaling  up  the  signal.  This  is  accomplished  by  solving  for  the 
desired  a  of  the  noise  (8) : 


SNR  =  10  LOG 


LIM 

£  sig“ (n) 

n=l _ 

LIM  2 

Z  noise“(n) 
n=l 


taking  the  antilog: 


(SNR/10) 


LIM  , 

Z  sig  (n) 

n=l _ _ 

LIM  2 

Z  noise  (n) 
n=l 


LIM 

Z  noise  (n) 
n=l 


LIM  ,f 
Z  sig“(n) 
n=l _ 

(SNK/10) 


LIM 

7 

Z  noise~(n) 

n=  1 _ 

LIM 


LIM  „ 

Z  sig  (n) 

n=l _ 

(SNR/10) 

LIM  •  10 


LIM 

9  1 

Z 

sig" (n) 

n=  1 

LIM  • 

,  „(SNR/10) 

10 

Therefore,  multiply  the  noise  samples  by  the  factor 


input 


_  _  # 
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FIGURE  2.  BLOCK  FLOW  DIAGRAM  OF  LATTICE  (INVERSE)  FILTER  IMPLEMENTATION 


T*mH-<F£*  FUNCTION 


Each  state  of  the  system  (see  Fig.  2)  for  an  impulse  function  input  to  the  one-formant  case. 


function  input  to  the  five-formant  case 
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RELATIONSHIPS  BETWEEN  LATTICE  PARAMETERS,  REE.  1  SYMBOLS,  AND  FORTRAN  VARIABLES 
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|  Case 

Pre  filter 

Prefilter 

Prefilter 

Lattice  Filter 

1  Figure  No. 

Formant,  Hz 

Bandwidth,  Hz 

y  k=  1 ,  ,  ,  .  Npo  1l* 

•\k,  k=l ,  .  .  .Nzere 

3 

2500.00 

150.00 

0.000000 

1 

0.000000 

0.985112 

0.995760  i 

4 

2500  .00 

150.00 

0.000000 

i 

-0.002642 

0.985112 

0.985765 

5 

650.30 

94.10 

-0.266366 

. 

-0.232250  ; 

1075.70 

91.40 

-0.522137 

-0.543236 

2463.10 

107.40 

0.174152 

0.14  5o  36 

3558.30 

198.70 

0.247624 

0.211085 

4631 .30 

89.80 

0.513600 

0.515390 

i 

0.260513 

0.260130 

0.141718 

0.067095 

-0.496971 

-0.420906 

-0.241914 

-0.221051 

0.943518 

0.860902 

6 

650.30 

94.10 

-0.266366 

-0.263710 

1075.70 

91.40 

-0.522137 

-0.529765 

2463.10 

107.40 

0.174152 

0.172208 

3558.30 

198.70 

0.247624 

0.254260 

4631 . 30 

89.80 

0.513600 

0.508698 

0.260513 

0.264767 

0.141718 

0.146636 

-0.496971 

-0.505169 

-0.241914 

-0.236822 

0.943518 

0.941123 

_ 

TABLE  II 

FORMANT  FREQUENCIES,  THEIR  BANDL'IDTHS,  AND  CORRESPONDING 
FILTER  COEFFICIENTS  FOR  THE  FOUR  INPUT  CASES  DEPICTED  IN  FIGS.  3-6 
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Command  Procedure  LATTICE.COM 

Controls  Sequence  of  Program  Execution  for  VT100  Terminals 


npii.v/usi  ii-i.»Pvs<:rsr.ii»/i,in.s»ssi.TPP»PYsr.r.«»Tioo/i.P> 


Command  Procedure  LATTEK.COM 

Controls  Sequence  of  Program  Execution  for  Tektronix  Terminal 


t  f*  / 1  i  *•«  ftp  \ i  *  /t.  r  n  ,  s  v<:<  c .  t  mb  t  u  v  •  itcstf  *  / 


28 


November  12,  1982 
BAC-LHS :cac 


jJ.«4H0  1s  4S‘llidn|  W  w  d  juo  d  4Ht  au  elili.lj  *.il  J  (Xtl'xiid.il  laiiU't  MlaS’lUjiji '*s)i  f  <*dU  i 

it  I'jiia*. 


Program  FILTER. FOR 

Recursive  Realization  of  Prefilter  and  Lattice  (Inverse)  Filter  Transfer  Functions 
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APPENDIX  B 

SAMPLE  TERMINAL  SESSION  FOR  THE  CASE  ILLUSTRATED  IN  FIG. 
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